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We extend our recent derivation of the time evolution equations for tiie energy content of magnetic 
fields and turbulent motions for incompressible, homogeneous, and isotropic turbulence to include 
the case of non- vanishing helicity. These equations are subsequently numerically integrated in order 
to predict the present day primordial magnetic field strength and correlation length, depending on 
its initial helicity and magnetic energy density. We find that all prior analytic predictions for helical 
magnetic fields, such as the epoch when they become maximally helical and their subsequent growth 
of correlation length L ~ a^'^ and decrease of magnetic field strength B ~ a~^'^ with scale factor a 
are well confirmed by the simulations. An initially fully helical primordial magnetic field is a factor 
4 X 10'* stronger at the present epoch then its non-helical counterpart when generated during the 
electroweak epoch. 
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I. INTRODUCTION 

The time evolution of primordial magnetic fields from 
very early times to the present is a subject of continuing 
research (for reviews see [l|-lJl)- It has become of par- 
ticular importance after the recent observational claim 
that a large fraction of the Universe may be filled by 
a magnetic field. A volume filling cosmic magnetic field 
is a natural prediction when magnetogenesis occurred in 
the early Universe. It is thus of interest to have as pre- 
cise as possible predictions on the evolution and possible 
magnitude of such fields. 

There has already been much effort put into this in- 
vestigation |6l-[l3l|. However, there is an important funda- 
mental difficulty which needs to be overcome, in particu- 
lar current numerical simulations do not reach the reso- 
lution required in order to reliably study magnetic fields 
over the large range of length scales and immense cosmic 
expansion factor between the early epoch of magnetoge- 
nesis and the recombination or present epoch. After all, 
the Universe changes in size by a factor of 10^^ between 
the electroweak and the recombination epoch while the 
respective dynamical time scales differ by a factor of 10^^. 
Only a small fraction of the initial magnetic energy sur- 
vives dissipation, how much exactly may though depend 
on the large scale tail of the magnetic field spectrum as 
has been seen to be the case for non-helical fields. Such 
tails arc extremely difficult to deduce from numerical sim- 
ulations. 

Though fairly detailed analytic estimates for the evo- 
lution of primordial magnetic fields exist |10l4l2l | , it may 



be of interest to verify them by more sophisticated meth- 
ods. One such method is a closure theory of magnetohy- 
drodynamics which has been recently applied by us |14{ 
to non-helical magnetic fields, with the important result 
that the field strength tail dependence on the scale L 
of B{L) ^ L~"/^ where a = 5 (cf. also [15[) seems to 
be essentially independent of the initial conditions in in- 
compressible MHD. Here the tail is important as it has 
a large influence over how much magnetic energy sur- 
vives. In this paper we extend the formalism developed 
in [IJI to include helicity. Helical ma gnet ic fields have 
already attracted much prior attention J15l - l2]| as helical 
magnetic fields may be much stronger due to an inverse 
cascade. 

To do so this paper is structured as follows: in Sec. [Tl] 
we give a short derivation of the time-evolution equa- 
tions for a homogeneous and isotropic medium including 
all terms which appear due to a non- vanishing magnetic 
helicity and apply the findings to the situation in the 
early Universe. More details can be found in Appendix 
[K\ In Sec. mil we then present our results of the time 
evolution of primordial magnetic fields according to our 
equations for different initial conditions and finally draw 
our conclusions in Sec. IVl 
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For an incompressible fiuid (i.e., V • v = 0) the equa- 
tions for the time evolution of the two main observables 
of magnetohydrodynamics, the velocity field of the tur- 
bulent medium v and the magnetic field B, are given 



by 



and 



9tB = - — AB + Vx(vxB) (1) 

Ancr 



av.-,v^V,v+<^^^ + f,„ ,2) 



and "Hq the magnetic helicity density, i.e 

d3 



hB^J^ A(x) . (V X A(x)) = I ^ A(x) . B(x) 

z p dkHk 



* ' d'fc <! ( p X B(k) ) • B(k)' 



(6) 



respectively. Here a is the conductivity (which is as- 
sumed to be very large in the following) and p the mass 

density of the fluid, while f^ is some viscous density force. where A is the magnetic vector potential and hs is the 
Note that for most part of the evolution of the early Uni- total magnetic helicity density, and where V • A = 0, i.e., 
verse incompressibility is an excellent assumption due to Coulomb gauge, has been assumed, 
the large speed of sound in the relativistic plasma 
Our main interest lies in the average build-up of mag- 



For all expressions V denotes the volume and we have 



1,.^. i-i i-1 n assumed cosmic homogeneity and isotropv which implies 

netic and kinetic energy density, respectively, as well as , ^^rr , ^ , r ■ i r i • i 

r.^ ,. , ,• •: , , 1,1 r ■ that Aife. (ife and Tit; arc lunctions only ot the magnitude 

of the magnetic helicity taken over an ensemble of cosmic , r , . ■, n,i n i,mi 



realizations (denoted by chevrons ()), i.e., 

{dtMq), (dtUq), {dtUq) 



(3) 



with Mq being the magnetic spectral energy density de- 
fined through 



k of the wave vector k. Furthermore, Parscval's Theo- 
rem has been used in all three cases in order to obtain 
a k integral where a hat denotes the Fourier Transform 
normalized by y 2 (cf. Appendix A of Ref. [iJI)- With 
these assumptions we find 



--^/t^^^W-/^!^^'^)!^^^/^^^^^- 



Uq the kinetic spectral energy density given by 



(4) 



T~La 



eK ^ f |^v2(x) = ^Jd^kmkr^pJdkU, 



M, = f^|B(q)p, 
;7, = 2V|v(q)P, 



qxB(q) •B(q)* 



(7) 
(8) 
(9) 



(5) 



By performing the calculations which are presented in Ref. [l4| and Appendix |X] of this paper, we obtain a very 
general result for the homogeneous and isotropic case which for ([3]) is given by 



(dtMq) 



I dfc|At< 



|g2 (M,) {Uu) - ^q' {Mq} {Mk) + Iji^l'k^ {Uq) {Hk) 



d9 



^ I4 (9' + fc' - qk cos 9) sin^ 9 (Mk) {Uk, ) 



(10) 



and 



(BtU,) = 



rH 
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Iq^k 



At{ - |g2 (Mfc) {Uq) - ^q^ {Uq) {Uk) + 



d9 



\ ^ {qk sin2 9 + 2kl cos 9) sin 9 {Mk) {Mk, ) 



31,2 



4 kl 



{3k~qcos9)sin-'9{Uk){Uk,) + 



1 q-^k 

(16^)2 kf 



-2q-q sin^ 9 + 2k cos 9) sin 9 {Hk) {Uk, ) 



as well as 

{dtUq) = 



r°° f ( 4 4 2 

J dkiAtI -e{Mq){nk) - -q^{Mk){Hq) - ^q'{Uk){nq) + 



d9 



^l'k\,.S 



2 kf 



sin' 9 {Uk,) {-Hk) 



(11) 



fl2) 



where q, k and fci are the magnitudes of the wave vec- tors q, k and ki = q — k, respectively, and 9 is the angle 




t3^' 


10-" 


&■ 






10-" 




10-2'' 




10-^5 



lO'* 10'' 

^/Mpc-' 



10' 



^/Mpc 



FIG. 1. Time evolution of magnetic (solid) and kinetic (dashed) spectral energies, as well as of the spectral magnetic helicity 
(dotted). Gray lines denote the initial conditions (i.e., at a = 1) while the black lines represent the situation for a = 10®. In 
the left panel maximal helicity has been assumed while in the right panel the case for vanishing initial helicity is shown for 
comparison. 



between q and k, i.e., q • k = q fc cos 9. Eqs. ([T0|) -([T2 |) are 
a set of well-defined equations since they ensure conser- 
vation of energy, momentum, mass, and helicity density 
to the lowest non-trivial order in At. 

It is important to note here that whereas the ideal 
MHD equations conserve energy and helicity, in any real 
application to astrophysics and cosmology dissipation 
plays an important role. In particular, the early Universe 
is described by a large Prandtl number |ll| , a measure of 
the relative importance of kinetic dissipation versus mag- 
netic diffusion. Dissipation is thus dominated by shear 
viscosity, such that magnetic diffusion can be neglected. 
Since energy dissipates by either kinetic or magnetic dis- 
sipation, but helicity only by magnetic dissipation, in the 
early Universe helicity is essentially conserved whereas 
energy is not. To model this in our simulations we intro- 
duce a viscosity term into Eq. (|lip . 

In addition it should be pointed out that the result 
obtained in |14l | is consistent with the equations presented 
here as it can be seen that for vanishing initial magnetic 
helicity it will be zero for all times and therefore simply 
the third term in ([TU)) and the last term in ([TT|) have to 
be dropped. 

In order to study the time evolution of magnetic fields 
in the early Universe we include expansion accounted by 
the scale factor a with og = 1 at the initial magneto- 
gensis epoch. In Section III of Rcf. [IJ], following also 
Ref. [nl, it was shown that Eqs. dTU]), (HI]), and P^ 
remain invariant in the expanding Universe when time 
derivatives d/dt are replaced by scale factor derivatives 
{Hq/q) 9/9 In a and physical quantities in the integrands 
of Eqs. (|10p - (|12p are replaced by their comoving ana- 
logues. These comoving quantities (marked by a 'c') are 
related to the physical ones by die = dta~^, Vc = w. 



Pc 



pa , 
Hqa' 



kc = fca, il/f^ — Mqa , U^ — Uqa , and 
^ whereas Hp, is the Hubble constant at the 



initial epoch. In the following the index 'c' is dropped 



and all quantities are meant to be the comoving ones 
unless noted otherwise. 



III. RESULTS 

It is now possible to analyze the time development of 
both the spectral magnetic and kinetic energy as well 
as the spectral helicity by numerically integrating the 
cosmic version of Eqs. (I10| - (|12p . In comparison to the 
results without helicity |14{ it may be seen from Fig. [T] 
that the absolute value and location of the integral scale 
(i.e., the scale at which the peak of the spectral energy 
is located) changes dramatically since due to an inverse 
cascade large amounts of energy are transferred to small 
q. On the other hand, the scaling of the large-scale tail 
with the wave number stays the same, i.e., we still have 
qMq ~ q^ and, in addition, we find qHq ^ q^ as well. 
In contrast to the non-helical case, however, the slope 
of the large-scale tail is much less important than in the 
maximally helical case the evolution of magnetic energy 
density and correlation length is determined solely by 
the requirement of helicity conservation, and not by the 
tail. In the following we present analytical arguments to 
explain these results. 



A. Large-Scale Magnetic Tail 

As in the case for non-helical fields [Ij] the non- 
vanishing term with lowest power in the Taylor expansion 
in {q/k) of Eq. (fT0|) - (fT2|) is of fourth order for large scales 
such that also for helical fields we have B ^ L^^ and 
V ~ L~2 on large scales. This result is also confirmed 
by simulation as can be seen in Fig. [TJ However, in con- 
trast to non- helical fields, in the fully helical field case, 
equipartition, i.e., Uk — Mk on all scales, does not hold 
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FIG. 2. Time evolution of the integral scale [upper left) and the values at the integral scale of the magnetic spectral energy 
(upper right), the kinetic spectral energy (lower left), and the total normalized magnetic helicity density hs/p = J Hkdk (lower 
right) for different values of the initial spectral helicity in fractions of its maximal value, Ho.max = SnMo/ko. Here we assumed 
that magnetogenesis has happened at the QCD epoch. 



anymore. Rather it is Uk — M^ only on small scales with 
Uk « Mk on the integral scale, an effect also observed 
inRef. [13. 



B. Evolution of the Integral Scale 

As it will be shown in the following, the main effect 
of non-vanishing initial helicity is the dramatic change 
of the time evolution of the integral scale and the corre- 
sponding values of the important quantities. 

First, we give a short summary of the time evolution 
of the magnetic spectral ener gy w ith vanishing initial he- 
licity. It has been shown [III, [23 that in this case, if the 
spectral energy has a slope proportional to fc"~^ (a > 1) 



for k < kj, i.e. 



{Mk 



Mo I — 
ko 



(13) 



where Mq is a normalization constant given by the value 
of Ml at ao = 1 , we have 



koa °+2 



(14) 



for the time dependent integral scale in the early ra- 
diation dominated Universe (i.e. H = Hq/u^), with 
fco = kj{ao) being determined by [22l | 



fco 



2wHo 






(15) 



where Uq — C//(ao). Since all the initial magnetic field 
has been dissipated away for k > ki(a), the magnetic 
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FIG. 3. Time evolution of the values of the magnetic and kinetic spectral energies and the helical spectral density at the 
integral scale, named Mi, Ui and Hi, respectively. For convenience these quantities have been multiplied by ki (for Mi and 
Ui) and by fc|/(87r) (for Hi). The plots show the situation for different initial values of the helicity "Ho- maximal helicity, i.e., 
Ho — 'Ho,max {upper left). Ho = W~^Ho,max {upper right). Ho = lO~'^Ho,max {lower right) and no helicity, i.e.. Ho = {lower 
left). Note that for the case of maximal helicity kJHi /{8tv) has approximately the same value as kiMi and is therefore not 
visible in the plot. Here we assumed that magnetogenesis has happened at the QCD epoch. 



field left over at scale factor a is simply given by Eq. 
inserted into Eq. ((T5)) . such that we find 



(16) 



If. however, the initial helicity is not vanishing, it can 
be seen from Fig. [2] that it only has an influence on the 
magnetic and kinetic spectral energies if at the integral 
scale it has a value which is close to the maximal spectral 
helicity value Hj^max given by 



n. 



Stt 



Mi 



(17) 



However, the magnetic spectral energy decays faster than 
the spectral helicity in the early Universe and therefore 
even by starting with a rather small initial helicity it is 
inevitable that after some time it obtains its largest pos- 



sible value given by (|17p . Therefore there are two succes- 
sive regimes in the time evolution, the first (referred to 
by '1' in the following) where the effect of helicity is neg- 
ligible and the second (referred to by '2') where helicity 
has its maximal value and dominates the time evolution. 
One can estimate the time it takes for some initial 
conditions to evolve to the point where a transition from 
regime 1 to regime 2 takes place. Here we will focus 
on the magnetic spectral energy as it is the important 
observable in our case. We are starting with some initial 
conditions where for qq = 1 we have 



kiMj = koMo, 



kj-Hi 

Sir 



Stt 



= fokoMo ., 



(18) 



such that /o = 1 corresponds to maximal initial helicity 
and /o = means no initial helicity at all, i.e., "Ho = 
foT~Lo.max ■ If we now assume that, as seen in Fig. [31 the 



time evolution for both has a power law dependence on 
the scale factor a, i.e., 



kjMi oc a"^i , kjUi oc a 



XI 



(19) 



where, as motivated above, ^i > Xi > 0, and we can es- 
timate the scale factor atr of the transition from regime 1 
to regime 2 by solving the condition fc/M/ = kj'Hi/{8n), 
i.e., 



koMoa-^'' = fokoMoa-'^' 
for a which gives 



o-tr — /o 



(20) 



(21) 



From this we can therefore derive a lower bound on /o, 
such that for all /o larger than this lower bound fmin the 
field has become maximally helical by recombination at 
scale factor arec 



fn 



(22) 



Here the epoch of recombination is relevant, rather than 
the present, since, as shown in Ref. |ll| . during matter 
domination after the epoch of recombination only negligi- 
ble further processing of the magnetic fields occurs, pro- 
vided the magnetic fields arc sufficiently strong. Taking 
the values which have been obtained, i.e., /ii = ^^^ = ^ 

for a = 5 and xi — 0.32, for Urec = 10* (implying mag- 
netogenesis during the QCD epoch at T ^ 30 MeV) we 
therefore get fmin — 7 x 10"^". This implies that even 
the smallest initial helicities ultimately have an effect on 
the time evolution of the fields. Furthermore, also the 
integral scale itself changes its dependence on the scale 
factor when entering the maximally helical regime 2, i.e., 
ki oc 0""^ in regime 1, where, according to (|T4|) . we have 
Ki = —To I niatching well the numerically found value, 
and fc/ oc a~'^'^ in regime 2 where we have numerically 
obtained K2 — 0.66 and fi2 — 0.67. 

The values K2 and ^2, as well as xij may be deduced 
unambiguously by the argument of helicity conservation 



LjBj ~ pkjHi 



'S,0 : 



(23) 



where Lj — 27: /kj is the integral scale and we have 
assumed that this scale carries the bulk of the helic- 
ity. From Eq. ([23]) and Eq. (|14|) one finds immediately 
Xi = ^^jT^ — 0.286 for a = 5, to be compared to the 
numerically found 0.3. The values of K2 and ^2 are ob- 
tained by considering the turbulent cascading of energy 
from the integral scale to the dissipation scale. This pro- 
cess will decrease Bi and will increase Lj. How fast this 
can happen is given by the requirement that an eddy 
turnover can take place on scale L/ at epoch with scale 
factor a |ll| . in particular 



SI 



p Lja 



2^2 



H ■ 



P_ 



(24) 



It needs to be stressed here that the above is a condi- 
tion to be formulated with physical quantities as opposed 
to comoving ones. Since we have dropped throughout 
the index 'c', i.e., all quantities in the above equations 
are comoving, additional powers of the scale factors a 



ph 



enter through Lj = Lja and p 



nP^^ 



p/a (note that 



Sf/. 



Bi/y/p). Using these equations we find 



L, 



Loa 



2/3 



Bi ~ Boa 



-1/3 



(25) 



for the evolution of comoving integral scale and magnetic 
field, or equivalently 



kr 



koa'^/^ kiMi ~ koMoa^^/^ 



(26) 



such that K2 — 2/3 and 112 = 2/3 which are in extremely 
good agreement with the numerical simulations. To sum- 
marize the time evolution we may write 



ki{a) 



kga '^^ , a < atr 
kna^r^^^a^^^, a > atr 



(27) 



with Ki = — ^ and K2 ~ 2/3 and 



kiMi{a) 



koMoa ^^ , a < atr 
/co^'/oafr ~'^^a^''^, a>atr 



(28) 



with /ii = 2-^^ and p2 ~ 2/3, where atr is given by 
Eq. (|2ip with pi ~ 2^^^ and xi ~ '^^p2 ^'^^ ^^^ initial 
conditions given by Eq. ([T8|) . 



C. Enhancement of the Energy Content and the 
Magnitude of the Magnetic Field 

With these arguments at hand one can now estimate 
the strength of the magnetic field which, compared to 
the case without magnetic helicity, is larger. In order to 
make a reasonable estimate it is important to take into 
account the time evolution of both the integral scale ki 
and the corresponding spectral magnetic energy Mi in 
the two different regimes introduced in Sec. IIIIBl given 
by (|27p and ([28]) . Using the formula for the magnetic 
field from [11] , 



B{L) = {STTkpMk)^ = (2fcMfc)= Bo ., 



(29) 



where Bq ~ (47rp)2 is the effective magnetic field for the 
Alfven velocity va — 1, i.e., the speed of light, for the 
optimistic case of initial equipartition between radiation 
and magnetic energies. This case would correspond to 
Bq — 3 X 10^^ G. If we now acknowledge the integral 
scale to be the coherence scale of the magnetic field, then 
for the magnetic field strength at Lj = 27r//s/ we get 



B{Li) ^ i2kiMi)^ Bo. 



(30) 
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FIG. 4. Present day magnetic field strength having survived dissipation during the evolution in the early Universe as a function 
of initial Alfven velocity va.o and helicity Ho — fo'Ho,max for magnetogenesis occurring during the QCD epoch {left panel) 
and electroweak epoch [right panel), respectively. Here all initial conditions above the dotted line lead to a maximally helical 
present day field. Note that es — pv\/2. It is stressed that the magnetic field has to fulfill va/L ~ H at the magnetogenesis 
epoch. For different initial conditions the modified va,o and /o to be used in the figure may be deduced from Eq. (f32ll . 



Using Eq. ()28|) we may then finally estimate the magni- 
tude of the present day magnetic field depending on its 
initial value and the helicity 



B, 



today 



3 X lO-^Gauss 



-5/7 „ ^ -8/7 
VA,0O,rec , JO S Clrec 

.1/3 -1/3 
VA.OJo '^rec , 



/o > a 



8/7 
rec 



(31) 



unless /o > 1 or VA,a > w^". Here w^% Lini, and fini 
quantify the initial total magnetic field strength, coher- 
ence length and helicity. Note that a = 5 should be 
employed in Eq. ((32|) as such a large scale magnetic field 
tail will quickly develop independent of the initial condi- 
tions 111. 



where va.o is the initial Alfven velocity at magnetogen- 
esis, i.e., t;^ Q ~ 2eB/ p, and where we have assumed the 
preferred exponent a = 5. The surviving magnetic field 
strength is enhanced in the helical case by a^ec /q over 
that in the non-helical case, which amounts to a factor 
2.7 X 10'^ /o and 3.7 x 10^/q for magnetogenesis during 
the QCD and electroweak epoch, respectively, highlight- 
ing the importance of helicity. The detailed prediction 
for the present magnetic field strength as a function of 
the initial magnetic field strength and helicity is shown in 
Fig.|3]where the QCD and electroweak epochs have been 
assumed as the epoch of magnetogenesis, respectively. In 
this figure it is assumed that va/L ~ H at the magne- 
togenesis epoch. If this is not the case, the appropriate 
initial conditions to be used in Fig. 2] may be obtained 



VA,0 



tm T (q+2) / ^A 
-'a ^mi 



H 



(32) 



/O — Jini \ „ 
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(71+21 



IV. CONCLUSIONS 

In this paper we have extended our recent deriva- 
tion [ij] of the evolution of the spectral magnetic and ki- 
netic energy densities in homogeneous and isotropic mag- 
netohydrodynamics by a closure theory to include helic- 
ity. The resulting ordinary differential equations were 
subsequently numerically integrated to predict the evo- 
lution of helical cosmic magnetic fields from the early 
Universe to the present. We find that the energy content 
and coherence length of sub-maximally helical magnetic 
fields essentially evolve as those of magnetic fields in the 
completely non-helical case. Due to the long cosmic ex- 
pansion from the early Universe to the present, however, 
even magnetic fields with the smallest initial helicities 
ultimately become maximally helical before the present 
epoch, thereafter decaying slower than non-helical mag- 
netic fields due to an inverse cascade. Our simulations 
agree well with all prior analytic estimates on the evolu- 
tion of helical magnetic fields. 
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Appendix A: Derivation of the Master Equations 
(flOl). pTj) and (fT^ 

The derivation of these equations follows exactly the 
same lines as the one presented in Appendix A of Ref. [1J| 
for non-helical magnetic fields. In the more general case 
presented here magnetic field correlators are assumed to 
be 



B„(k,i)Bb(k',t')*) ^ C7i5kk"5«'A4. Sa. 






+ C3d]ik'Stt'ieabckc'Hk 



(Al) 



with the first and second part of the correlator describing the non-helical and helical components of the magnetic 
field, respectively. 



Appendix B: Alternative form of the Master Equations 

We give here an alternative formulation of the master equations (|10p - (|12p in terms of ki = q — k which is more 
suitable for numerical integration: 



rq+k 






At< 



r^q^ (M,) {Uk) - ^q' {Mq) {Mk) + Ij^l'k' (H,) {Hk) 



dfci 



q — k\ 



q 



q^k 3q^ 



9'fci 



qk^ 



qk qki qkf 



and 



mk^kf i6fc3fci wkkf mkf sfcfci i6fc3 mkf mki mk wk^ 



(Mk) (L/feJ 



(Bl) 



(dtUg) = J dkUtl- ^q' {Mk) {U,) - ^<z^ ([/,) (C/fe) 

q^ q^k 3g^ qk^ 3qk bqki 



q+k 



dfci 



q-k\ 



16kkf 8kf 8kki Wkf 



Iq^ 



3q^ 



ci 16fc 

.3 



{Mk){Mk,] 



llq^k 5q^ Sq^ki 5qk^ 9qk iqki qkf \ 
32k^kf 32kkf 32k^ki '' 32kf ^ 16fcfci ^ 32k^ ^ 32kf ^ 32ki ^ 32k ^ 32 k^ ) 



(B2) 



(Uk) (C/feJ 



I f q^ 



3q^k q^ki bqk^ 3qkki qkf 



(87r)2 \16kki 
as well as 

{dtnq 



16/ci 



J dklAtl^k''{M,){nk 



^Kmxmj 



^q^AIk}{H,}-^q^Uk){n,} 



q+k 



dfci 



\q-k\ 



q^k 



qk^ qk qki 



8kkf ikf Akki 8kf 4fci 



{Uk,){Hk) 



(B3) 



We note here that for some terms a substantially shorter expression than that given in Appendix B of Ref. [1^ has 



been found by realizing that J^ dfc/|^_^,| dfci is equivalent to J^ dfci J, 



rq+k 
l\q-k\ 



dk. 
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